Macroscopic Superpositions of Phase States with Bose-Einstein Condensates 
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Quantum superpositions of macroscopically distinguishable states having distinct phases can be 
created with a Bose-Einstein condensate trapped in a periodic potential. The experimental signature 
is contained in the phase distribution of the interference patterns obtained after releasing the traps. 
Moreover, in the double well case, this distribution exhibits a dramatic dependence on the parity of 
the total number of atoms. We finally show that, for single well occupations up to a few hundred 
atoms, the macroscopic quantum superposition can be robust enough against decoherence to be 
experimentally revealable within current technology. 

PACS numbers: 67.85.-d 03.75.Gg 67.85.Hj 



Introduction. Nonlinearity is crucial for the creation of a 
superposition of macroscopically distinguishable states, 
often referred to in the literature as "Schrodinger's cat" 
[0, 0, 0, H|. It was first suggested by Yurke and Stoler 
|l| to use a kerr medium to create a superposition of 
photonic coherent states having different phases and 
to detect it via homodyne interferometry. However, it 
was quickly realized that the experimental realization of 
these special states is rather prohibitive due to the typi- 
cally weak optical nonlinearities and strong losses in the 
medium. Recent advances in atomic physics, where large 
nonlinearities arise naturally, are changing this scenario 
and encouraging a renewal of early proposals. Macro- 
scopic superpositions with a small number of trapped 
ions have been created in Q , while Bose-Einstein conden- 
sates (BECs) might provide even brighter sources thanks 
to their inter-particle interaction which introduce a large 
kerr-like coupling of matter waves. In particular, several 
thermodynamical and dynamical schemes for creating a 
superposition of two BEC states differing by a macro- 
scopically large number of particles have been proposed 
in the literature jf|. These so called "NOON states" max- 
imize particle entanglement and can be useful in quantum 
information protocols as, for instance, Heisenberg-limited 
inter ferometric phase estimation. 

In this manuscript we discuss an experimentally feasi- 
ble protocol for the creation and detection of a macro- 
scopic superposition of states having different relative 
phases with a Bose-Einstein condensate trapped in a pe- 
riodic potential. These entangled states [l| are generated 
by the nonlinear unitary evolution governed by the de- 
coupled TV^-mode Bose-Hubbard Hamiltonian: Hn s = 
J2f=i e jCjCj + ^^CjCjCj and are closely related, as will 
be shown later, to NOON states. Here c- (ct) annihi- 
late (create) a particle in the j th condensate, ej is the 
energy offset due, for instance, to an external potential 
superimposed to the optical lattice and Uj is the single- 
condensate interaction energy. For a double well system, 
Ns = 2, we show, by simulating the formation of sev- 
eral single-shot interference density profiles with a many- 
body Montecarlo technique [7[, that the relative phase 
probability distribution contains distinct peaks. Each 



peak corresponds to a phase state component of the su- 
perposition and its position depends dramatically on the 
parity of the total number of atoms. We demonstrate 
that a clear signature of the creation of these macro- 
scopic superpositions also appears in the interference of 
an array of BECs, even in a single-shot density profile. 
We finally account for the problem of decoherence due to 
one, two, and three-body losses [8[. For typical experi- 
mental trapping parameters and single well occupations 
up to a few hundred atoms, the decoherence time can be 
about 500 ms, longer than the typical formation time of 
the macroscopic superposition of phase states which can 
therefore be experimentally created and detected within 
current technology. 

Double well potential. We first consider a BEC trapped 
in a symmetric double well potential (ei = 62, IA\ = 
U2 = U) and rewrite the Hamiltonian in the pseudo 
angular momentum representation, H2 = UJ^ @, 0]. 
As initial condition we choose a state of Nt particles 
IV^o)> = En=o C ne- ineo \N T -n,n) z (where \N T -n,n) z 
is a number Fock state [10]). The state = 
e -iH 2 tn/ 2 /h^Q Q }^ obtained after a time = , can 
be formally written as 
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where £ = (£ = 1) for N T even (odd) [11]. It is in- 
structive to project the state Eq.([T]) over the SU(2) ba- 
sis states \N T ,0) = V2W V r+1 E^o^I^t - n,n) z 
[T^ . If the initial |^(^o)) has a relative phase distri- 
bution localized about #o, the phase distribution of the 
evolved state, P(0) = \(Nt, 0\^(t 7V / 2 )) | 2 ? is characterized 
by 2 peaks, see Fig. [TJ which is the consequence of be- 
ing 1^(^/2)) a superposition of 2 states having different 
relative phases. Since the peaks have a finite width (see 
discussion below), the highest visibility is reached just 
at when the peaks are maximally separated. For 

instance, when the initial state is given by a binomial 

distribution, C n = 2jv ^ /2 \J ( N j) and Oq = 0, the two 
components of the superposition are exactly orthogonal 
(^(f 0I^(-tt + f 0) = 0, and P(0) has peaks of width 
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w ~ 1 / yj Nt • With trapped BEC, a realistic scheme 
for the creation of a superposition of two states hav- 
ing different relative phases involves the sudden split- 
ting of a single condensate [Hj]. The BEC is left in a 
state slightly squeezed in the relative number of particles, 
hfo) = |^(0o = 0)) ~ EnZo e ~ (n ~ NT/2)2/ ^\ N T-n,n) z , 
which provides the initial condition of the decoupled non- 
linear evolution. The width of the relative number distri- 
bution is <j s — y/Nr/(2s) and s is the squeezing param- 
eter. After a time this state evolves in the super- 
position Eq.([T]) (for moderate initial number-squeezing 
the overlap between the two components is exponentially 
small). We release the confining potential and let the 
condensate ballistically expand and overlap, giving rise 
to an interference pattern [3, EH from which we ex- 
tract a single value of the relative phase. We will show 
that the phase distribution obtained upon several inter- 
ference experiments is reasonably well described by the 
SU(2) probability P(0), when the initial state is number- 
squeezed s > 1. Let us consider a simple model which 
can be numerically studied with a Montecarlo technique 
0- First, we simulate the formation of a single-shot in- 
terference density profile. The probability to detect Nt 
particles at the same time and positions {x± : #2, %n t } 
is 

(^ft f )|0t(x 1 )...0t(^ T )0( XjVT )..^( Xl )|^ft f )) 
{x} ~ N T \ 

Here <j>(x) = <j)\{x)c\ + ^2(^)^2, where <j)j{x) is a nor- 
malized wave function of the j th condensate. The nu- 
merical calculation is simplified by using two counter- 
propagating plane waves: <j>i{x) = 02 (%)* — e t7TX , with 
x G [0, 1] . After generating Nt random positions dis- 
tributed with P{^}, the phase is extracted by fitting the 
density profile to the function p(x; 0) = 1 + cos(27nr + 6). 
We repeat several times the interference protocol to ob- 
tain a probability phase distribution. In FIG. 0Ja),(b) 
we show the results of 400 independent phase estimations 
(dots) for Nt = 10 (a) and Nt = 11 (b). The agreement 
with the SU(2) phase distribution (solid line) holds al- 
ready for small values of Nt and improves for a larger 
number of atoms as long as s > 1. Thus, the single-shot 
density profile can be approximated by interfering spatial 
wave functions with a relative phase randomly sampled 
from P(9). We emphasize that this agreement is not ob- 
vious. In particular, SU(2) fails for phase squeezed states, 
s < 1, whose phase distributions are better matched by 
projecting over binomial phase states. 
In FIG. [2(a), (b) we can distinguish two peaks separated 
by 7r, each corresponding to a different phase component 
of Eq. (pQ) . The phase shift due to a change in the parity 
of Nt can be clearly seen as a shift of tt/2 among the 
distributions of (a) and (b). As expected, the width of 
the peaks w ~ s/^JNt increases with the relative num- 
ber squeezing of the initial state. In typical experimental 
conditions there is no control on the parity of the total 
number of atoms. In this case, the system is described 
by a classical mixture of superposition states, half corre- 
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FIG. 1: (color online). a),b) Montecarlo phase distribution 
(dots, with the phase values distributed over 20 bins) and the 
SU(2) phase distribution P{6) (solid line), for the state Eq.(pQ) 
with s = 1.5. In (a) Nt is even and P(0) ~ cos 2Nt (0); in (b) 
Nt is odd and P(9) ~ sin 2iVT (#). c),d) Schematic pseudo 
angular momentum representation of Eq.(pQ) for even (c) and 
odd (d) values of Nt- 



sponding to odd Nt and half to even Nt- For instance, 
the phase distribution at would show 4 peaks, with 
widths scaling as ~ 1/ ^Nt-, thus clearly distinguishable 
already for small Nt- Notice that the occurrence of two 
peaks (or four if the parity of Nt is not controlled) in 
the phase distribution constitutes, by itself, a signature 
of the presence of a quantum superposition, rather than 
a statistical (non coherent) mixture with the same com- 
ponents. Indeed, if the system had decayed into a mix- 
ture at any time t < t^li-* it could not have evolved into a 
macroscopic superposition. This because the state Eq.([T|) 
is created and lives in a very narrow temporal window, 
~t 7r/2 ±h/UN T 0. 

To conclude, we point out that, when the initial state is 
given by a binomial distribution, we can rewrite Eq.([T]) as 
|tf(t f )> = (e^A^O), + e i *\0,N T ) x )/V2 ( |tt(* f )> = 
(e-** \Nt, fyy+e^ |0, N T ) y )/V2 ) for Nt even (odd) 0. 
The representation of these states in the Bloch-sphere 

is shown in FIG. [2(c), (d). The even ( odd ) l*(*f)> are 
simply linked by a tt/2 collective rotation about the z 
axis, and may be both transformed into a NOON state 
(e-^\N T ,0) z + e^\0^T}z)/V2 with a tt/2 rotation 
about the y (x) axis [17|, [l8[. Therefore, a macroscopic 
superposition of two binomial states carries the same 
amount of particle entanglement (which is conserve d by 
local unitary transformations) of the NOON state [191 ]. 
As a consequence, it would be, for instance, as useful as 
the NOON state for quantum information protocols. We 
notice that, in the same way, a NOON state can be cre- 
ated with high fidelity from a slightly number squeezed 
initial state. 



Lattice potential The protocol for creating phase cats in 
the double- well has a direct generalization to an array of 
condensates. We consider a superfluid BEC trapped in a 
one-dimensional optical lattice modulated by a harmonic 
potential. The initial state, a product of coherent or 
slightly squeezed states Yl^i \ s j(®)) (sol localized in each 
lattice well, is created by increasing the interwell barriers 
rapidly enough to prevent the system from reaching the 
Mott insulating phase 0. The state then evolves with 
the decoupled Hamiltonian Hn s : 

l*i(*)>=C;5> ^ e-^-e-'^nr^ln), (3) 

where Nj is the mean occupation number of the j th site, 
Cj is a normalization constant, a s = ^/Nj / 's <C Nt, 
and s is the squeezing parameter which depends on 
the lattice ramping time [21]. At time t n the lattice 
is switched off and an image of the cloud is taken af- 
ter a time of flight tf. To construct the single-shot 
interference density profile p{x\ 6q, . . . , 0n s -i), we sam- 
ple the values of the phase of the wave-function in each 
well with the distribution Pj(0) oc \(0\sj(t n ))\ 2 , being 
1^ = "71^ 5^n5) e m6> |n). A simple calculation yelds 

AT , „ im( X -jd) 2 (x-jd) 2 

p(x;O u ...,e Ns ) <x lEfii^e^'e 2ht ' 
where d is the lattice period, m is the atomic mass, 
I = htf/mlo and Iq is the initial width of the Gaus- 
sian wave function describing each well along the lat- 
tice direction. Finally, we extract a single value for 
the phase by fitting the density profile to the function 
[1+/3 cos(# + 2ttx/ L)]G(x), where ft and are fitting pa- 
rameters, L = htf/md is the period of the first-harmonic 
modulation of the profile, and G(x) is a Gaussian enve- 
lope that accounts for the finite dimensions of the system 

FIG.[2Ka),(b) show typical single-shot spatial density pro- 
files of Nt = 1000 87 Rb atoms loaded into N s = 164 
sites of a ID optical lattice and squeezing parameter 
s = 3/2. FIG. [2ja) corresponds to a fitted value = 
while FIG. [2](b) to = tt. In FIG. E^c) we plot the 
polar diagram after 400 shots. Dots are fitting results 
of repeated simulations: the distance from the origin is 
the amplitude ft, while the polar angle is the phase 0. 
In FIG.[2Kd) we show the phase probability distribution: 
the values of = and = tt are strongly favored, thus 
revealing the presence of a macroscopic superposition of 
phase states. The loss of visibility, compared to the dou- 
ble well case, is caused by zero-point energy differences 
among neighboring sites as will be discussed below. In 
the lattice case faster oscillations with wavelength L/n 
( n integer ) can also be fitted out (23[. The latter cor- 
respond to the combination of relative phases between 
wells distant nd. Thus, for instance, phases fitted from 
two oscillations having distinct wavelengths, can differ 
by 7r (each fit can provide or tt with equal probability) . 
This would be a clear signature of the creation of the su- 
perposition obtained in a single interference experiment. 
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FIG. 2: (color online). a),b) Single-shot density profiles (blue 
solid line) for = and = tt, respectively. The yellow 
solid line is the fitting function (see text), c) polar plot and 
d) probability phase distribution (dots) after 400 interference 
shots. The blue solid line in d) is a guide to the eye and the 
phase values are distributed among 20 bins. Here Nt = 1000 
87 Rb atoms released from a ID optical lattice, with a time 
of flight tf = 20 ms. The trapping parameters are a = 5.3 
nm, d 0.27 um, V /h 243 KHz, cu± = 2tt x 70 Hz and 
uj x = 27T x 6 Hz, and the formation time of the superposition 
is t n ~ 60 ms. 



Discussion. The superposition of phase states of a few 
atoms have probably been already created experimen- 
tally during the collapse and revival of a matter wave field 
trapped in a three-dimensional optical lattice [24]. Un- 
fortunately, no signature could have been seen in the in- 
terference patterns, due to column averaging [22[. A one- 
dimensional configuration is therefore optimal. More- 
over, in order to avoid dephasing and drifting, the trap- 
ping configuration must be such that Uj and e^, respec- 
tively, take site independent values. In typical experi- 
ments the one-dimensional lattice is superimposed to a 
confining harmonic potential with frequencies uo x along 
the lattice direction and oo± in the transverse directions, 
which gives ej = Q(j — iVs/2) 2 , with Q = muld 2 /2. 
Therefore, we require the maximum drift rate between 
neighboring sites to be much smaller than the creation 
time of the phase states superposition, QNs/h <C l/t n . 
In order to have a site-independent interaction energy, 
we need the each local chemical potential fij = NjUj <C 
w x ,w±, where uo x is the frequency of each well along the 
lattice direction. It is worth to emphasize that these con- 
straints are much relaxed for a double well setup where 
any asymmetry between the wells would just lead to a 
global drift of the relative phase distribution, making the 
creation protocol very robust, while an unbalanced occu- 
pation of the wells would only decrease the visibility of 
the spatial interference fringes. 

Decoherence. Finally, we take into account the pos- 



sibility for one, two and three-body losses which 
can rapidly destroy the coherent superposition of atomic 
states. The condition to make these processes negligi- 
ble is \[ p H c < 1 (p = 1,2,3), where t c is the formation 
time of the macroscopic superposition of phase states, 
\[p] = K^Nj J drcj)f(r) is the loss rate relative to the 
p-body processand is an experimentally measured 
rate constant [25]. The optimal configuration is obtained 
for a particle density low enough to have only a few losses 
events (usually the three-body losses dominate since they 
scale as density to the third power), but high enough to 
allow a clear extraction of the phase from the interference 
pattern. For the lattice setup used in FIG. [2j we have 
X [ q ] U = 0.002, and \ [ * ] U = 0.0005, in the central well, 
with tn ~ 60 ms. On the other hand, for the experimen- 
tal double well trap parameters of [15[, for instance, each 
well should contain no more than Nj ~ 400 particles 
which would give a formation time approximately 
equal to the decoherence time, ~ 500 ms. 
Conclusions. A macroscopic superposition of phase 
states can be realized experimentally with a condensate 
trapped in a periodic potential. The macroscopic coher- 
ence is robust against asymmetries and decoherence and 
can be unambiguous detected from the interference pat- 
terns of the overlapping condensates. 
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